Vortex Viscosity in Magnetic Superconductors Due to Radiation of Spin Waves 
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In type-II superconductors that contain a lattice of magnetic moments, vortices polarize the 
magnetic system inducing additional contributions to the vortex mass, vortex viscosity and vortex- 
vortex interaction. Extra magnetic viscosity is caused by radiation of spin-waves by a moving vortex. 
Like in the case of Cherenkov radiation, this effect has a characteristic threshold behavior and the 
resulting vortex viscosity may be comparable to the well-known Bardeen-Stephen contribution. 
The threshold behavior leads to an anomaly in the current-voltage characteristics, and a drop in 
dissipation for a current interval that is determined by the magnetic excitation spectrum. 



Magnetic ordering coexists with superconductivity in 
many different compounds that are known as magnetic 
superconductors. For instance, the RN12B2C borocar- 
bides (T] or Ba(Fe 1 _ :c Co :! ;)2As2 (pnictides) [2] , are super- 
conductors that comprise a lattice of magnetic moments 
in their structure. Coexistence of magnetism and super- 
conductivity is also realized in hybrid systems that are 
specifically designed for this purpose [3] . A lot of work is 
being devoted to understanding the interaction between 
vortices and magnetic domain walls [3], with the main 
goal of enhancing vortex pinning. However, much less is 
known about the effect of the magnetic degrees of free- 
dom in the dynamic regimes of flux creep and flux flow. 

Vortices are topological defects of the superconduct- 
ing order parameter that emerge in type II superconduc- 
tors under the application of a strong enough magnetic 
field. In magnetic superconductors, the magnetic mo- 
ments couple to the magnetic field of the vortex via Zee- 
man interaction. Therefore, vortex motion in a magnetic 
superconductor is accompanied by a dynamic readjust- 
ment of the magnetic moments. Since the low-energy 
excitations of a magnetically ordered system are spin- 
waves or magnons, a moving vortex can emit magnons 
under certain conditions. Such magnon emission, in turn, 
affects the vortex motion. 

In this letter we demonstrate that the interaction be- 
tween vortices and magnetic moments of a type II mag- 
netic superconductor leads to a contribution to the vortex 
mass, vortex viscosity and vortex-vortex interaction. We 
show that a moving vortex can transfer energy to the 
magnetic moments by emitting spin-waves. This energy 
transfer gives rise to a "magnetic viscosity" . The inertial 
part of the response of the magnetic system modifies the 
dependence of the vortex energy on its velocity. Thus, 
it generates a magnetic contribution to the vortex mass. 
The viscosity of this " magnetically dressed" vortex may 
also differ significantly from the "bare" vortex value. The 
spin- waves have a well defined dispersion relation Sl(fc). 
Similar to the physics of Cherenkov radiation of photons 
by a moving electron, a superconducting vortex moving 
with velocity v radiates spin-waves with momentum k 
when the kinematic condition £l(k) = v • k is satisfied. 



As a consequence, the magnetic viscosity of a vortex in 
an antiferromagnet (or a ferromagnet with easy-axis per- 
pendicular to the applied magnetic field) has a threshold 
behavior in the vortex velocity, i.e., it is effective only if 
v exceeds a critical velocity v c . 

The radiation of photons in the classical Cherenkov 
problem takes place when a charged particle moves faster 
than the speed of light in a dielectric medium. In a 
magnetic superconductor, the threshold velocity v c de- 
pends on the spin-waves velocity, the size of the vortex 
core, and the gap of the spin-wave excitation spectrum. 
The threshold behavior of the magnetic viscosity, rj m (v), 
can be observed experimentally as an anomaly in the I-V 
characteristics of the flux flow regime, near the voltage 
that corresponds to the threshold velocity v c . An im- 
portant result is that r) m (v > v c ) may be comparable in 
magnitude to the core viscosity, rj ci which is induced by 
excitations of the quasiparticles in the vortex core [7]. 
Here we only consider temperatures that are well below 
the magnetic transition temperature T/v- 

The physics of the magnetic mass is analogous to the 
polaron problem. [ 9 The inertia of the magnetic moments 
adds up to the vortex inertia. In this paper we are mainly 
concerned with the vortex flow regime for which the ef- 
fects of the vortex mass are not important. [TT] 

1. Spin-vortex interaction in a magnetic superconduc- 
tor. The magnetic field of a vortex j, oriented along the 
z axis and located at the spatial coordinate Rj in the xy 
plane, has the form zH*(r — Rj), where z is the unit vec- 
tor along the z-axis. The Fourier components of H*(r) 
are H* k = (f> /(l + \ 2 k 2 ), where A is the London penetra- 
tion length and <po is the flux quantum. In a quasi-static 
approach (y much smaller than the Fermi velocity), the 
vortex moving along the trajectory Rj(t) induces a time- 
dependent magnetic field if* [r — Rj (t)] that acts on the 
magnetic subsystem. Application of a uniform magnetic 
field induces a uniform internal field Hq = zHq that 
creates a vortex lattice. The magnetic field inside the 
superconductor, H = ziP(r), has the periodicity of the 
vortex lattice (in an isotropic superconductor the vortices 
are oriented along the direction of the internal magnetic 
field.) A moving vortex lattice creates a time-dependent 
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magnetic field H(r) = Hy[r ~ Rj'M]- For a dense 
vortex lattice, Hq (j)o/X 2 , the periodic component of 
the magnetic field is small compared to the uniform com- 
ponent, h z (r) = H z (r) - H < H . 
The Zeeman term is: 

H int = -J2 ^ Z [v n - RMSnit), (1) 

nj 

where r„ are the coordinates of the spins described by 
the operators S n and \x — gfis (<? is the gyromagnetic 
factor and \ib is the Bohr magneton). We describe spins 
quantum mechanically, while the vortices are treated as 
classical objects. The vortex-spin interaction given by 
Eq. (JTJ leads to the following equations of motion for 
spins and vortices 

ihS n = [S„, H M ] fih z (r n - Rj)[S„, S*] - ihTS n , 

j 

(2) 

rj c Rj - --j~2^^ [oj +b h + hp 

+ ^ F v (RjRj/), (3) 
j' 

where [A, B] = AB — BA and L z is the system length 
along the z axis, T is the spin relaxation rate due to 
coupling of spins with other degrees of freedom, and ry c is 
the vortex Bardeen-Stephen core viscosity. Fl and F p are 
the Lorentz and pinning forces respectively, F v (Rj, Rj') 
is the non-magnetic force acting on vortex j due to the 
presence of other vortices located at R 3 <. These forces 
are defined per unit of vortex length. (...) indicates the 
quantum mechanical average. The magnetic system is 
described by the Hamiltonian, Hm, which is specified 
below and accounts for exchange coupling between spins, 
anisotropy, and Zeeman coupling to Ho. 

The first term in the right hand side of Eq. ^ describes 
the interaction of a given vortex with other vortices due 
to polarization of the spin system. When vortices are at 
rest, this term may lead to a change in the vortex lattice 
structure depending on the characteristics of the mag- 
netic subsystem and on the applied magnetic field. For 
dynamical properties, this term leads to additional dissi- 
pation due to energy transfer to spin excitations, vortex 
mass renormalization caused by the inertia of the spins, 
and change of the moving vortex structure. In the fol- 
lowing we only consider dynamical effects of the perfect 
square vortex lattice, Rj(i) = R° + vt, with reciprocal 
lattice vectors G = (2tt /b)(m x , m y ), where m x y are inte- 
gers, n v = 1/b 2 is the vortex density, and b it the vortex 
lattice parameter. We assume that spins are close to 
equilibrium, i.e., that T is strong enough in comparison 
with the rate of spin excitations. The periodic compo- 
nent of the magnetic field h z (r, t) is weak in comparison 
with Hq. Thus we use the linear response approxima- 
tion with respect to h z (r,t) in order to find the Fourier 



components (S^(lu)). We note that only nonzero spa- 
tial harmonics, G / 0, contribute to the coupling of the 
moving vortex lattice to the magnetic subsystem. The 
spatial Fourier components of h z (r,t) at the initial time 
t = are h z G {0) = (/) n v /[l + A 2 G 2 ]. We obtain: 

n(Sh<t)) = -^g(O) f W*(G,i-t')e-* G - vt ',(4) 
n v Jo 

z n 

Here a is the lattice parameter of the magnetic ions 
and \ zz (G,t) is the longitudinal differential susceptibil- 
ity of the spin system with respect to the alternating field 
h G {t). By expressing (S z (uj)} in Eq. ^ via Eq. Q, we 
obtain the equation of motion for the vortex lattice: 

7? C R= — V |^(0)| 2 Ge lG ' R ( t) f dt' X zz (G,t- t')e- lG - vt ' 

n v G JO 

+ F L , (6) 

where R is the center of mass coordinate of the vortex 
lattice. Here we neglected the pinning force. When the 
vortex lattice moves with constant velocity v, we obtain 
from Eq. ^ : rj c V + F m = F L and 

F m = — ^ E I&g(0)| 2 G ■ v hn.\ X "(G,u = G • v)]. 

" G 

(7) 

Here we introduced the magnetic viscous force F m (per 
unit of vortex length). So far the discussion is general 
and valid for any magnetic superconductor. 

Magnetic viscous force for an antiferromagnet. We 
now consider the simple case of a two-sublattice anti- 
ferromagnet with an easy axis anisotropy along the x 
direction. Thus, we introduce the following magnetic 
Hamiltonian for spin S moments: 

(i,n) (l,n) n 

that is defined on a bipartite lattice with coordination 
number Z. Here (I, n) indicates that I and n are nearest- 
neighbor (NN) lattice sites, J is the exchange constant 
between NN, and Ja > is the amplitude of the ex- 
change anisotropy term. This model is a good starting 
point for describing magnetic superconductors that have 
a Neel temperature, Tjy, lower than the superconducting 
critical temperature T c . While the exchange term sta- 
bilizes a two-sublattice antiferromagnetic ordering along 
the x direction (easy axis), Ho induces a uniform com- 
ponent along the z-axis. This canted antiferromagnetic 
phase persists at T -C Tjv until the field reaches the sat- 
uration value H s , i.e., the staggered magnetization com- 
ponent (J,(S*) — ±[iSy/l — H 2 /H 2 vanishes for H > H s 
(different signs correspond to different sublattices). 



3 



We will compute the longitudinal susceptibility, 
% zz (k, lu), in the spin- wave approximation. Since the 
natural upper cut-off for the momenta of the vortex lat- 
tice field is l/£ <C 1/a, we compute x zz (k,w) in the 
long wave-length limit. The result is x zz (k, ui) = for 
H > H, and 



X zz (u,,k;# ) = 



Vt 2 + s 2 k 2 



«2 



(w + iry 



(9) 



. The saturation field is H s = (fi/a 3 )S/xo, 
X zz (u = 0, k = 0) is the static suscepti- 



(v/ha 3 / 2 )s^jXz(T 



Hi/Hi) and 



for H < H s 
where xo = 
bility. Here f2i 
si = (n/haV^Sy/Jil - H$/Hf). The dispersion tl k of 
the spin-wave excitations that couple to the vortex field 

is [g 
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Here fi and so are the spin wave gap and velocity. For 
H > H s , the staggered magnetization is zero and the 
radiation of spin- waves vanishes at low temperatures be- 
cause the magnetic field cannot change the already sat- 
urated amplitude of the uniform magnetization. 
By inserting Eq. ^ into Eq. §f§, we obtain: 



2n v v 2 
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G • v r(o^ + s\G 2 ) 



n G [(G 
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r 2 ]' 



(ii) 

To proceed, we replace the sum over G by an integral 
by assuming that T/s >• 27r/&, i.e., that the integrand 
varies over a length scale much longer than 2ir/b. When 
r <C k*v, with k* = VLq/ yjv 2 — Sq, the Lorentzian can 
be replaced by a delta function, Im[x zz (k, to — k ■ v)] ~ 
(il 2 + s 2 k 2 )6(flk — k • v)/fife. In this form the kinematic 
condition of the radiation is explicit. As we mentioned 
above, linear response is applicable when the spin relax- 
ation is fast compared to the rate at which magnons are 
created by the vortex lattice. We can now specify this 
condition explicitly: T^>a 3 / 2 fli(j) /(2iTyn^X) 2 . By per- 
forming the angular integration, k • v = kv cos <f>, we ob- 
tain the magnetic viscous force per unit of vortex length 
for v > sq: 



F m (v)-- 



A 2 



L kdk nl + sjk 2 9{k-k*) 



o (2^) 2 y^fc 2 - n 2 [k 2 + A 



-212 



,(12) 



where 6(x) = 1 if x > and zero otherwise. Unlike the 
Cherenkov electron, which is a point-like object and gen- 
erates a field with arbitrary high wave vector, the vortex 
has a finite core size and its spatial Fourier components 
are exponentially suppressed for k > 27r/£. This deter- 
mines the upper integration limit in Eq. (|12[). 



We will now analyse this equation for F m (v) in the 
range: 



27tsq a 3 / 2 fii / r/> ( 



/On 



(2ttA) : 



2irso 



(13) 



As the vortex velocity increases, the value of k*(v) starts 
at k*(so) — oo, for v = sq, arL d decreases monotonically 
for v > sq. Depending on the value of k*(v) relative to 
the upper integration limits, Eq. ( 12 ) has two different 
regimes: 

(i) k* is above the upper integration limit, k* > 2tt/£. 
The magnetic viscous force vanishes, F m — 0, for v < 
v c =[s 2 + (& /2ir) 2 ] x / 2 . 



(ii) For k* < 2n/^ the integration in Eq. (12 1 has the 
limits k* < k < 2n /£. At v — v c the dissipation function 
has a threshold behavior F m oc y/v — v c . Note that the 
threshold is sharp because T <C 2ttsq/^. As the velocity 
increases beyond [v 2 + (^ilo/27r) 2 ] 1 / 2 , but remains small 
enough to ensure k*(v) ^> 2ir/b, the viscous force has an 
approximately linear dependence on the vortex velocity 
and we can define a magnetic viscosity via 
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Vn 



■(v 2 -s 2 + 2s 2 /xo)- 



Xo ( 0o 
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?V 

(14) 

Fig.JTJshows the dependence of the magnetic viscous force 
F m on the vortex velocity. As H approaches H s , the am- 
plitude of the periodic part of the magnetic field decreases 
and the value of F max eventually vanishes. 

As temperature increases and approaches T/v, the spin- 
wave relaxation rate T increases. The dynamics becomes 
relaxational near the magnetic transition[8l I15j. In this 
regime, the magnetic viscosity does not exhibit a thresh- 
old behavior in the vortex velocity. 

3. Experimental signatures. We can estimate the mag- 



nitude of the magnetic viscosity by evaluating Eq. (14 1. 
For typical values of the relevant parameters, fi ~ 10/is, 
T c ~ 10 K, H c2 ~ 10 T, p n ~ 100 ^ficm, A - 200 cm, s ~ 
50m/s and Ho = 2 K, we obtain rj m ~ 10 -7 g/(cm • sec). 
This value is comparable to the BS core viscosity [7] 
Vc = (-ffc20o cr n)/c 2 ~ 10~ 7 g/(cm • sec). Clearly, the 
magnetic contribution to the vortex viscosity is more im- 
portant for nearly isotropic magnetic materials (small 
gap Qq) such as superconducting compounds or het- 
erostructures whose magnetic moments are provided by 
transition metal ions. 

We will now discuss the current- voltage characteristics. 
In the experimental set up, the voltage V is measured at 
a fixed current. The electric current exerts a Lorentz 
force, Fl = (<j>o/c)I, on the vortex. The voltage V is 
proportional to the velocity of the vortex flow via the 
Josephson relation V = (h/2e)(v/b). The vortex veloc- 
ity is determined from the balance between the Lorentz 
force, Fl, and the viscous force given by Eq.([7|. As the 
current exceeds the value I OIL at which the vortex veloc- 
ity (determined by the Bardeen-Stephen viscosity only) 
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FIG. 1: (Color online) (a) Total viscous force in arbitrary 
units as a function of v/so, as given by Eq. (12 1. The blue 



dashed line represents the effect of Bardeen-Stephen viscosity 
due to vortex core alone, (b) Effect of the magnetic viscous 
force on the I-V characteristics. Magnetic viscous force results 
in a voltage drop relative to the BS result for currents above 
7 on corresponding to the vortex velocity v c . 



equals to the onset velocity, v c , the character of TV curve 
changes. Right above v c (onset of region II in Fig.[l]), the 
magnetic contribution to the vortex viscosity is nonzero 
and increases fast, F m oc ^Jv — v c . Consequently, vortex 
velocity is weakly dependent on / near the threshold and 
V - V c oc (I - I on ) 2 where V c = {h/2e)(v c /b). Once the 
current is large enough to propel vortex to a larger speed 
where magnetic viscosity saturates to a constant value, 
the linear behavior, V cx I, is recovered. The slope of 
linear behavior of V vs I in region II is smaller compared 
with region I because the total viscosity is larger. 

It is important to determine if the critical vortex ve- 
locity v c , necessary to reach the onset of magnetic dis- 
sipation, is experimentally attainable. The limiting fac- 
tors are heating and the Larkin-Ovchinnikov instability 
of the vortex flow at high velocities. [Ml [17] A current 
density of ~ 50 KA/cm is needed to reach vortex veloc- 
ities close to the spin-wave velocity, s ~ 50m/s. If heat- 
ing is the limiting factor, such currents can be reached 
with pulsed measurements. In addition, s decreases for 
increasing magnetic field and temperature. It is hard 
to get a reliable estimate of the critical velocity for the 
Larkin-Ovchinnikov instability, which is given approxi- 
mately by v* vp^T e \/T- m (vp is the Fermi velocity and 
Tel.in are the elastic and inelastic scattering times respec- 
tively). Measurements of YBCO films yield v* ~ 1000 



m/s [T7]. 

In conclusion, moving vortices radiate spin-waves in 
superconducting antiferromagnets for v > v c . This effect 
decreases the flux flow resistance, i.e., it reduces energy 
dissipation. The onset of magnetic viscosity appears in 
the TV characteristic curve as a local deviation at the 
threshold voltage that is determined by a vortex veloc- 
ity v c . For materials with T/v < T c , we predict a drop 
in the resistance as temperature drops below T/v. The 
same effect can be observed in hybrid structures made 
of alternating superconducting and magnetic layers for 
a magnetic field perpendicular to the layers. The mag- 
netic layers can be antiferromagnetic or ferromagnetic 
with easy axis parallel to the layers. 
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